Shear banding and interfacial instability in planar Poiseuille flow 
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Motivated by the need for a theoretical study in a planar geometry that can easily be implemented 
experimentally, we study the pressure driven Poiseuille flow of a shear banding fluid. After discussing 
the "basic states" predicted by a one dimensional calculation that assumes a flat interface between 
the bands, we proceed to demonstrate such an interface to be unstable with respect to the growth 
of undulations along it. We give results for the growth rate and wavevector of the most unstable 
mode that grows initially, as well as for the ultimate flow patterns to which the instability leads. We 
discuss the relevance of our predictions to the present state of the experimental literature concerning 
interfacial instabilities of shear banded flows, in both conventional rheometers and microfluidic 
channels. 

PACS numbers: 47.50.+d, 47.20.-k, 36.20.-r. 
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I. INTRODUCTION 

Complex fluids have internal mesoscopic structure that 
is readily reorganised by an imposed shear flow. This 
reorganisation in turn feeds back on the flow field, re- 
sulting in strongly nonlinear constitutive properties. In 
some systems this nonlinearity is so pronounced that 
the underlying constitutive curve relating shear stress 
Txy to shear rate 7 in homogeneous flow is predicted 
to have a region of negative slope dT^y/d'y < [U, 
In this regime, an initially homogeneous flow is unsta- 
ble to the formation of coexisting shear bands of differ- 
ing local viscosities and internal structuring, with band 
normals in the flow-gradient direction y [3]. The signa- 
ture of this transition in bulk rheometry is the presence 
of characteristic kinks, plateaus and non-monotonicities 
in the composite flow curve Q. Explicit observation 
of the bands is made using local rheological techniques 
such as flow birefringence [5] and NMR d, 0], ultra- 
sound d, , heterodyne dynamic light scattering 
or particle image llj velocimetry. Using these meth- 
ods, the existence of shear banding has been firmly es- 
tablished in a wide range of complex fluids, including 
n,JIl E ES [13, El, El and 
25| surfactants; side-chai n liq - 
291; viral suspensions [13, l28| |: 
polymer 



wormlike 
lamellar [2; 

uid crystalline polymers 

telechelic polymers [29]; soft glasses [30, 31^^ 
solutions [33] ; and colloidal suspensions Q . 

Beyond the basic observation of shear banding, ex- 
periments with enhanced spatial and temporal resolu- 
tion have more recently revealed the presence of complex 
spatio-temporal patterns and dynamics in many shear 
banded flows AM^MM^ [H [M, 31, M M, \M, [M [12, 




49[. In many such cases, the bulk 
stress response of the system to a steady imposed shear 
rate (or vice versa) is intrinsically unsteady, showing ei- 
ther temporal oscillations or erratic fluctuations about 
the average (flow curve) value. Local rheological mea- 



surements reveal such signals commonly to be associated 
with a cornplicated behaviour of the interface between 

the bands g M, M M, M, M\M,\M,\M,M,\M,\M- 

The majority of these measurements have been in one 
spatial dimension (ID), normal to the interface between 
the bands. However 2D observations in Refs. [13, \^ 
explicitly revealed the presence of undulations along the 
interface, in a boundary driven curved Couette flow, ac- 
companied by Taylor-like vortices [4§]. The undulations 
were shown to be either static or dynamic, according to 
the imposed flow parameters. 

Theoretically, instability of an initially flat interface 
between shear bands was predicted in boundary driven 
planar Couette flow in Refs. 0, 151], |5^. In this 
work, separate 2D studies in the flow/flow- gradient {x- 
y) Ho, [EH and flow-gradient/ vorticity {y-z) [s^] planes 
revealed instability with respect to undulations along the 
interface with wavevector in the flow and vorticity direc- 
tions respectively. In both cases the mechanism for insta- 
bility was suggested to be a jump in normal stress across 
the interface 53]. 

While these predictions provide a good starting point, 
there remains the possibility that the interfacial undula- 
tions observed in Refs. [13, [H, [4^ originate instead in 
curvature driven effects such as a bulk viscoelastic insta- 
bility of the Taylor Couette [54| type in the high shear 
band, as discussed in Ref. [4^. These were neglected in 
the planar calculations of Refs. [13, [5l|, [13] (Other pos- 
sibilities, also neglected, include free surface instabilities 
at the open rhcometer edges; and an erratic stick-slip 
motion at the solid walls of the flow cell. We shall not 
consider these further in what follows here either.) 

In principle, therefore, either (or both) of (at least) two 
possible mechanisms could underlie the observed interfa- 
cial undulations: (i) a bulk viscoelastic Taylor Couette 
like instability of the strongly sheared band, or (ii) in- 
stability of the interface between the bands, driven by 
the normal stress jump across it. Of these, scenario (i) 
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FIG. 1: Left) Dotted line: homogeneous constitutive curve 
for a = 0.3, 77 = 0.05. Thick sohd hne: composite flow curve 
for one-dimensional planar shear banded Couette flow (data 
already published in Ref. [s^l). Selected stress Tsei = 0.506. 
Thin solid line: parametric plot of local shear stress T^yiy) — 
^xy{y) + vyiy) = Cd — y) against local shear rate 7(2/) for 
one- dimensional planar shear banded Poiseuille flow with I = 
0.00125, G = 2.0. 

Right) Halved pressure gradient versus total throughput for 
one- dimensional planar Poiseuille flow with a — 0.3, rj = 0.05, 
I = 0.0025. This shows a kink at the onset of banding at 
G/2 = Tsci as expected. 

can only arise in a curved geometry. Experiments in pla- 
nar shear could therefore in principle help discriminate 
between these scenarios, by eliminating the curvature re- 
quired for (ii). However they are technically difficult to 
implement in a boundary driven setup. 

There thus exists a clear need for theoretical predic- 
tions in a planar flow geometry that could easily be im- 
plemented experimentally. An obvious candidate com- 
prises pressure driven flow in a rectilinear microchan- 
nel of rectangular cross section with a high aspect ratio 
Lz/Ly ^ 1. Indeed, such experiments have recently been 
performed [5l,[5l,[53l, as discussed in more detail below. 
With this motivation in mind, in this paper we study 
the planar Poiseuille flow of a shear banding fluid driven 
along the main flow direction a; by a constant pressure 
drop dxP = —G. For simplicity we assume the fluid to 
be sandwiched between stationary infinite parallel plates 
at y = {0, Lj,}, neglecting the lateral walls in the z direc- 
tion, and so taking the limit L^jLy — > cx) at the outset. 
Our main contribution will be to show an interface be- 
tween shear bands to be unstable in this pressure driven 
geometry, as it is in the boundary driven planar Couette 
flow studied previously [13, HH, [12] . We will furthermore 
give results for the growth rate and wavevector associ- 
ated with the early stage kinetics of this instability, as 



well as the ultimate flow patterns to which it leads. 

The paper is structured as follows. After introducing 
the rheological model and boundary conditions in Sec. [Til 
we calculate in Sec. IIIII the one-dimensional (ID) shear 
banded states that are predicted when spatial variations 
are permitted only in the flow gradient direction y, arti- 
ficially assuming translational invariance in x and z, and 
accordingly assuming a flat interface between the bands. 
These form the "basic states" and initial conditions to be 
used in the stability calculations of the rest of the paper. 

In Sec. IIVI we study the linear stability of these ID 
basic states with respect to small amplitude perturba- 
tions with wavevector qj.x in the fiow direction. As in 
the case of boundary driven flow studied previously, we 
find an undulatory instability of the interface between 
the bands [13, [EH, HI] . Results are then presented for the 
ultimate nonlinear dynamical attractor in this x-y plane, 
from simulations that adopt periodic boundaries in x. 
This exhibits interfacial undulations of finite amplitude 
that convect along the flow direction at a constant speed. 
In Sec. |V] we turn instead to the flow-gradient/vorticity 
plane y-z, likewise demonstrating linear instability of the 
interface with respect to small amplitude perturbations 
with wavevector q^i. We also give results for the ulti- 
mate nonlinear flow state, which in this plane is steady. 
Directions for future work, which will include full 3D cal- 
culations, are discussed in Sec. IVII 



II. MODEL AND GEOMETRY 

The generalised Navier-Stokes equation for a viscoelas- 
tic material in a Newtonian solvent of viscosity 77 and 
density p is 

p{dt + v.V)v = V.(T - PI) = V.(S + 2^T) - PI), (1) 

where v(r, t) is the velocity field and D is the symmet- 
ric part of the velocity gradient tensor, (Vv)^^ = davp. 
Throughout we will assume zero Reynolds' number p = 0. 
The pressure field P{r,t) is determined by enforcing in- 
compressibility, 

V-v = 0. (2) 

The quantity S(r, t) in Eqn. [T] is the extra stress con- 
tributed to the total stress T(r,t) by the viscoelas- 
tic component. We assume this to obey the diffusive 
Johnson-Segalman (D JS) model [sl, [5§| 



(a* -f V • V) S = a(D • S + S ■ D) + (E • O + O • S) + 2GoD - 



-V^S. 



(3) 



r 



Here a is a slip parameter, Gq is a plateau modulus, r metric part of the velocity gradient tensor. The diffusive 
is the viscoelastic relaxation time, and Jl is the antisym- 
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term V^S is needed to correctly capture the structure 
of the interface between the shear bands, with a sUghtly 
diffusive interfacial thickness 0{l), and to ensure unique 
selection of the shear stress at which banding occurs [60] ■ 
Within this model we study flow between infinite flat 
parallel plates at y = {0, Lj,}. The fluid is driven in 
the positive x direction by a constant pressure gradient 
dxP — —G, the plates being held stationary. At the 
plates we assume conditions of zero flux normal to the 
wall n.VSa^ = V a, /? for the viscoelastic stress (al- 
though other choices are possible [HH), with no slip and 
no permeation for the fluid velocity. Throughout we use 
units in which Go — 1,t = 1 and Ly = 1. Unless other- 
wise stated we use a (dimensionless) value rj = 0.05, sug- 
gested by the experiments of Refs. [59]. We take a = 0.3, 
although our results are quite robust with respect to vari- 
ations in this quantity. An order of magnitude estimate 
suggests I = 0(10-3) [el. 
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FIG. 2: Left) Velocity profiles in one-dimensional planar 
shear banded Poiseuille fiow for a = 0.3, r] = 0.05, I — 0.0025. 
Increased throughputs correspond to increased applied pres- 
sure gradients G = 2.0, 3.0, 4.0, 5.0, 6.0, 7.0, 8.0, 9.0, 10.0. 
Right) Corresponding shear rate profiles, which follow as the 
spatial derivative of the velocity profiles. 



III. ONE-DIMENSIONAL BASIC STATE 

In this section we discuss the flow curves and shear 
banded states predicted by ID calculations that allow 
spatial variations only in the flow-gradient direction y, 
assuming (often artificially) translational invariance in 
the flow direction x and vorticity direction z. As a warm- 
up discussion to the case of planar Poiseuille flow (PPF) 
that forms the primary interest of this paper, we review 
first the more familiar case of planar Couette flow (PCF). 

In steady ID PCF, the shear stress T^y is uniform 
across the flow cell. This follows trivially from solving 
Eqn. [1] in this geometry. Within the (often artificial) 
assumption of a similarly homogeneous shear rate field, 
corresponding to a velocity field v = jyx, the consti- 
tutive relation is then given by Txy{j) = T,xy{'j) + rjj, 
where T,xy{'y) follows as the xy component of the solu- 
tion of Eqn. [3] obtained by assuming stationarity in time 
and homogeneity in space. See the dotted line in Fig. [T] 
(left). For an applied shear rate in the region of nega- 
tive slope dTxy/dj < 0, homogeneous flow is predicted 
to be linearly unstable with respect to the formation of 
shear bands. The steady state composite flow curve is 
then shown by the thick solid line in the same figure. 
For shear stresses T^y < Tsei (resp. T^y > Tgd), where 
the selected stress Tsc\ = 0.506 for the particular choice 
of parameters in Fig. [U the system shows homogeneous 
flow on the low shear (resp. high shear) branch of the 
constitutive curve. Applying a value of the shear rate 7 
in the window between these branches, we obtain shear 
bands that coexist at the selected stress Tj;y — Tg^i- 

In steady ID PPF driven by a pressure drop dxP = 
—G, Eqn. [1] dictates the shear stress to be linear across 
the flow cell: T^y = —G{y — i). Well away from any 
interfaces (the almost vertical regions in Fig.[2l right) the 
value of the shear rate 7 then follows as the lower (resp. 
upper) solution branch of the same relation Txy{'j) = 
^xy{'j)+il'i discussed above for PCF, in the regions where 



\Txy\ < Tsci (resp. \Txy \ > Tgci). This leads to the steady 
flow states shown in Fig. [21 which are shear banded for 
all values of the applied pressure drop G > 2Tsci shown. 
In bulk rheology, the signature of the transition to shear 
banded flow is the pronounced kink seen in the composite 
flow curve of Fig. [T] (right). 

The ID shear banded flows discussed in this section 
will form the basic states and initial conditions for the 
stability studies of the rest of the paper. In Sec. IIVI 
we will consider 2D flow in the x-y plane, with periodic 
boundaries in the flow direction x, and assuming transla- 
tional invariance in z. In Sec. [Vl we will turn to 2D flow 
in the y-z plane, with periodic boundaries in the vorticity 
direction z, and assuming translational invariance in x. 



IV. FLOW, FLOW-GRADIENT PLANE 

In this section we relax the assumption of transla- 
tional invariance in the flow direction x and perform 
a two-dimensional study in the flow/flow-gradient {x- 
y) plane. As before we consider flow between closed 
walls aX y = {0, Ly}, driven by a constant pressure drop 
dxP — —G. The boundaries in x are taken to be periodic. 
For simplicity, translational invariance is still assumed in 
the vorticity direction z. 

Our study comprises two parts. In Sec. IIV Al we study 
the linear stability of the ID shear banded basic states of 
Fig. [2] with respect to perturbations that have infinitesi- 
mal amplitude and a wavevector ^ajX in the flow direction. 
We will show these basic states to be linearly unstable to 
such perturbations in most regimes. An initial condi- 
tion comprising shear banded flow with a flat interface is 
thereby predicted to evolve towards a 2D state that has 
modulations along the interface, via the growth of these 
perturbations. In Sec. IIVBI we study the model's full 
nonlinear dynamics in this x-y plane, giving results for 
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FIG. 3: Top) Dispersion relation of growth rate lu versus 
wavevector Qx during the initial stage of 2D instability, start- 
ing with a ID shear banded planar Poiseuille basic state. 
In each case a = 0.3, -q — 0.05, G = 4.0. Circles: ex- 
tracted from early time dynamics of the full nonlinear code. 
Groups of data upward correspond to I = 0.02, 0.015, 0.01, 
0.005. Within each group, data are shown for {Dt, Nx, Ny) — 
(0.000025,200,400), (0.0000125, 200, 400), (0.000025, 400, 
800). These are mostly indistinguishable, demonstrating 
convergence with respect to grid and timestep. Dashed 
lines; cubic splines through the data for {Dt,Nx,Ny) = 
(0.000025, 200,400), as a guide to the eye. Lx = 4.0 in each 
case. Solid line: results of analytical calculation in the true 
limit I — > 0. Crosses: extrapolation of the I ^ results to 
I — using the scaling uj{qx, I) = u){qx,l = 0) — a{qx)l. 
Bottom) Growth rate (left) and wavevector (right) at the 
peak of dispersion relations as a function of the pressure drop 
G, calculated analytically in the limit / 0. 



the ultimate 2D flow state that is attained once these 
modulations have grown to, and (as we shall demon- 
strate) saturated at, a finite amplitude. 

Before presenting our results we discuss briefiy our nu- 
merical method. In previous work the full nonlinear dy- 
namics of the DJS model in boundary driven PCF were 
simulated by SMF in the flow/flow-gradient plane [63j . 
assuming translational invariance in the vorticity direc- 
tion. During the course of that study, the code was 
carefully checked against an earlier calculation of the 
linear stability of an interface between shear bands in 



FIG. 4: Greyscale snapshots of the order parameter 'Exxix, y) 
at a representative time t = 60.0 on the ultimate nonlinear 
attractor for a = 0.3, rj = 0.05, I = 0.005 for Lx = 4.0 and 
G = 2.0, 4.0, 6.0, 8.0 (top four subfigures downwards) with 
Dt = 0.000025, Nx = 200, Ny = 400. Bottom subfigure is for 
G = 4.0 with a longer cell Lx = 6.0 (so Nx = 300). 

PCF [13, HH, providing a stringent check that the non- 
linear code was working correctly. A straightforward two- 
line modification adapts that code to the pressure driven 
PPF of interest here. 



A. Linear stability analysis 

In each run of this modified code we start with an ini- 
tial condition comprising a ID banded state of Sec. IIIIl 
corresponding to shear banded flow with a flat interface 
between the bands. To this state we add 2D pertur- 
bations of tiny amplitude. By monitoring, in this fully 
nonlinear code, the early-time growth of the Fourier com- 
ponents e'K'p{iqxx) exp(a;t) of the perturbation, we can 
extract the dispersion relation uj{qx) that characterises 
the linear instability of the ID basic state. These nu- 
merical results for this dispersion relation are shown by 
dashed lines in Fig. [3] (top), for different values of the 
interfacial thickness I. As can be seen, at any fixed value 
of qx the growth rate shows a linear dependence on I: 
uj*{l,qx) — = 0,qx) — a{qx)l where the intercept 

uo*[l = Q,qx) plotted versus qx accordingly forms the 
dispersion relation extrapolated to Z = 0, shown by the 
crosses in Fig. [3] (top). The value / = cannot be ac- 
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FIG. 5: Top) Snapshot of interface height h{x) for a = 0.3, 
rj = 0.05, I = 0.005, = 4.0 and G = 2.0 at six times 
equally spaced by At = 0.05265. Dt = 0.000025, A^^ = 200, 
Ny — 400. Inset: the same data plotted versus transformed 
coordinate x — ct revealing simple convective motion with a 
constant wavespeed c. 

Bottom) Diamonds: wavespeed c versus pressure drop G, 
obtained by performing the transformation shown in Fig. [5] 
a = 0.3, 77 = 0.05, I = 0.005, = 4.0 for runs with Dt = 
0.000025, = 200, Ny = 400. Circles: fluid velocity at 
interface, taken from the ID banded proflles of Fig. [2] 



cessed in this full nonlinear code, because this limit is 
pathological in shear banding fluids. However, analytical 
linear stability calculations can nonetheless still be per- 
formed in this limit |5]| . The solid line in Fig. [3] (top) 
accordingly shows the results of a truly linear analytical 
stability calculation, performed by HJW in the true limit 
I ^ Q. As can be seen, this indeed agrees well with the 
extrapolation of the numerical results to / = 0. Because 
the I ^ and I = calculations were performed indepen- 
dently by the two authors, and using different methods, 
this provides a good cross check between both sets of 
results. 

Because of the linear nature of this stability system, 
any eigenvector can be separated into sinuous and vari- 
cose contributions, in which the perturbation to the flow 
field across the channel is odd (snake-like) and even 
(sausage-like) respectively. Each of these is an eigen- 
vector in its own right. In the I = calculations, this 
separation is carried out explicitly. We find that the sin- 



uous (snakelike) modes are always more unstable than 
varicose ones. Because our interest is in the most un- 
stable mode, the / = results presented here are all for 
sinuous perturbations. (In our runs of the full code for 
I ^ Q, we do not impose any such symmetry a priori. 
Instead, it emerges naturally from the full system.) 

The location of the peak in the dispersion relation is 
shown as a function of pressure drop G in Fig. [3] (bot- 
tom), for the / = calculation. As can be seen, for 
values of G well within the banding regime the growth 
rate uj = 0(1), comparable to the linear viscoelastic re- 
laxation time, and the wavelength A^; = lixjqx — 0(2), 
comparable to twice the gap width. This is consistent 
with the presence of secondary velocity rolls [i^, |4§| of a 
wavelength comparable to the gap width, and also con- 
sistent with the dominant wavelength seen in our non- 
linear results discussed in Fig. |4] below. (Such velocity 
rolls will be shown explicitly in our corresponding study 
of the flow-gradient/ vorticity plane in Fig. [7] below.) As 
the applied shear rate approaches the low shear branch 
from above, and the width of the high shear bands ac- 
cordingly tends to zero, the growth rate and wavelength 
of the most unstable perturbations also tend to zero. 

For a small range of pressure drops 1.20 < G < 1.52, 
two peaks are evident in the dispersion relation (not 
shown). For G < 1.3 the peak with the larger value of qx 
is the more unstable, with crossover to dominance of the 
lower Qz peak for G > 1.32. This causes a "kink" in the 
plot of a;* against G at G « 1.3, only just discernible in 
Fig. El 

B. Ultimate nonlinear dynamics 

At long times the instability described above saturates 
at the level of finite amplitude undulations along the in- 
terface. Representative greyscale snapshots of the sys- 
tem's eventual state are shown in Fig. |4l As can seen in 
Fig. [5] (top) , the interfacial undulations convert along the 
positive X direction with constant speed. This speed is 
plotted as a function of pressure drop in Fig. [5] (bottom), 
and is comparable, but not exactly equal, to the value of 
the fluid velocity at the interface. As noted above, our 
I code does not impose any sinuous/ varicose symme- 
try. Indeed, the ultimate nonlinear attractor can be seen 
in Fig. [4] to be dominated by a mode that is neither fully 
sinuous not fully varicose. 



V. FLOW-GRADIENT, VORTICITY PLANE 

In this section we turn attention to the flow- 
gradient /vorticity (y-z) plane, now assuming transla- 
tional invariance in the flow direction x. As usual we 
assume closed walls in y, and (for simplicity) periodic 
boundaries in z. We will return in the conclusion to dis- 
cuss briefly the implications of having focused only on 
2D studies in this paper. 
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FIG. 6: Top) Dispersion relation of growth rate u versus 
wavevector during the initial stage of 2D instability, start- 
ing with a ID shear banded planar Poiseuille basic state. In 
each case a — 0.3, -q — 0.05, G = 4.0. Thin solid lines: calcu- 
lated by linear stability analysis using the method of Ref . [121 , 
with lines upward corresponding to I — 0.0035, 0.003, 0.0025, 
0.002, 0.0015. Thick solid lines: extrapolation of the Z / 
results shown by each of the thin solid lines to Z = using the 
scaling io{qx, I) = (-o{qx, 1 = 0)— a{qx)l- Indistinguishability of 
the thick lines demonstrates the validity of this scaling. Cir- 
cles: extracted from early time dynamics of the full nonlinear 
code for / = 0.0035. 

Bottom) Growth rate (left) and wavevector (right) at the 
peak of dispersion relation as a function of the pressure drop 
G for Z = 0.003, 0.0025, 0.002, 0.0015 (sets upwards in left 
hand figure). 



A. Linear stability analysis 

We consider first the linear stability properties of the 
ID basic states of Fig. [5] with respect to fluctuations with 
wavevector in the vorticity direction, exjp(iqzZ + ujt). To 
do so we perform a truly linear calculation, calculating 
the eigenmodes of the stability matrix that is obtained 
by linearising the full nonlinear equations about the basic 
state. This linearisation problem was in fact studied pre- 
viously in the flow-gradient/ vorticity plane in the context 
of planar Couette flow in Ref. [H^l • The same linearised 
code can be used here, simply inserting the new pres- 
sure driven basic state as an input to the elements of the 
matrix. 



The resulting dispersion relations u}{qz) are shown in 
Fig. E] (top). As is evident, at each value of q^ the growth 
rate increases linearly with decreasing values of the inter- 
facial width Z, allowing us to extrapolate to the case I = 
(thick line in Fig.[6l top), as in the case of the x-y plane 
above. The location lo* , k* of the peak in the dispersion 
relations is plotted as a function of the pressure drop G 
in Fig. [51 bottom. As can be seen, for large values of the 
interfacial width I the instability is absent (all oj* < 0). 
For intermediate values of I there is a window in G be- 
tween the onset of shear banding and the onset of this 
interfacial instability. The size of this window decreases 
with decreasing Z, and we believe would extrapolate to 
zero as Z ^ 0, such that all shear banded states are lin- 
early unstable in the limit of a thin interface. 

These truly linear results were checked for consistency 
against the linearised dynamics of a fully nonlinear code, 
which directly simulates the full 2D dynamics of the 
model in the yz plane. (The main use of this code is 
to generate the results in the next subsection.) In each 
run of this nonlinear code we used as an initial condi- 
tion the ID basic state of Fig. [2 subject to small 2D 
fluctuations in the y-z plane. By monitoring the early- 
time growth of the Fourier components exp{iqzz) exp(a-'t) 
we extracted the dispersion relation oj{qz) for the linear 
(in)stability of the ID basic state. As shown by the circles 
in Fig. [51 good agreement is found with the true linear 
stability calculation described in the previous paragraph. 
This provides a stringent cross check between our linear 
stability calculation and our nonlinear code. 



B. Ultimate nonlinear attractor 

In previous work we simulated the full nonlinear dy- 
namics of the DJS model in 2D plane Couette flow in the 
flow-gradient/ vorticity {y-z) plane assuming trans- 
lational invariance in the flow direction x. A straightfor- 
ward two-line modification allows this code to be adapted 
to the case of pressure driven flow of interest here. As de- 
scribed in the previous section, we performed a series of 
runs of this code starting with a ID basic state, subject 
to small 2D perturbations. We consider in this section 
the ultimate state that the system attains at long times. 
In each case, the instability was found to saturate in a 
steady state with finite amplitude undulations along the 
interface (Fig. [71 top). Such undulations have recently 
been imaged experimentally in the pressure driven flow 
of wormlike micelles in microchannels of high aspect ra- 
57| . Associated with these undulations is a secondary 
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flow comprising the velocity rolls shown in the second 
subflgure (Fig. [71 bottom). 



VI. SUMMARY AND OUTLOOK 

In this work, we have studied the pressure driven pla- 
nar Poiseuille flow of a shear banding fluid sandwiched 
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FIG. 7: Order parameter 'Ea:xiy,z) grayscale (top) and veloc- 
ity map showing secondary flow in y-z plane (bottom): snap- 
shot at time t — 1895 in steady state for a — 0.3, rj = 0.05, 
I = 0.0015 for L:, = 4.0 and G = 2.0 with Dt = 0.05, 
= 1000, Ny = 1066. 

between infinite stationary flat parallel plates. Using a 
combination of linear stability analysis and direct numer- 
ical simulation, we have shown an initially flat interface 
between shear bands to be unstable with respect to the 
growth of undulations along it. The early time growth 
rate of these undulations scales as the reciprocal stress 
relaxation time of the fluid, and the corresponding wave- 
length is comparable to the separation of the plates. At 
long times the undulations saturate in a finite amplitude, 
cutoff by the nonlinear effects of shear. 

Throughout our study we have neglected any in-ffow 
and out-flow effects at the start and end of the channel, 
assuming a well established flow field that is free from 
end effects. Associated with this is our assumption that 
the spatial limits relevant to this established fiow field 
correspond to the temporal limit t — > co in each run 
of our code. We have furthermore neglected the effect 
of any lateral walls in the z direction, taking from the 
outset the limit Lz/ Ly oo. In practice one might ex- 
pect lateral walls to impose a further loss of translational 
symmetry in the z direction, as explored in Ref. [STj . 
Perhaps most importantly, in conducting separate two- 



dimensional studies in the x — y and y — z planes, we 
have neglected the possibility that exists in three spatial 
dimensions of nonlinear interactions between the Qx and 
qz modes. Indeed, although the Qx modes have much 
faster growth rates than the qz modes in the early time 
(linear) growth regime, they are cut off more aggressively 
at nonlinear order by the effects of shear. This results in 
comparable amplitudes for the qx and qz models in the 
ultimate attractors of our separate 2D studies. It there- 
fore seems unlikely that 3D effects cannot be neglected, 
as we shall explore in future. We mention finally that 
the constitutive model used in the present work is highly 
oversimplified, in particular in having a Newtonian high 
shear branch. 

Recent experiments in microchannel flow used ID ve- 
locimetry imaging to reconstruct the local flow curves of a 
shear banding wormlike micellar surfactant solution [s^ . 
Surprisingly, these curves were found to be non-universal, 
failing to collapse when plotted for different applied pres- 
sure gradients on a single plot. Possible explanations of 
this observation, to be explored in future work, include: a 
larger ratio 1/ Ly than considered here; coupling between 
flow and concentration; and the effects of secondary 3D 
flows. 
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